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We discuss the recent progress in computing the D-meson decay constant and D — > ir£v form 
factors from the lattice QCD simulations with Nf = 2 dynamical using Wilson quarks. We 
report f D = 201(20) MeV and F+(l GeV 2 )// D = 4.04(78) GeV -1 at a ~ 0.08 fm. 

1 Introduction and lattice setup 

Accurate determinations of the Cabibbo-Kobayashi-Maskawa (CKM) couplings provide an es- 
sential test of the Standard Model. The most straightforward method for their extraction is 
through the leptonic and/or semileptonic meson decays. However, such an extraction from the 
experimently measured decay widths requires a reliable information about the hadronic quan- 
tities, namely the decay constants and the form factors. Such an information is expected to be 
provided by the QCD simulations on the lattice. 

Many lattice calculations of these quantities have been attempted over the past twenty 
years. Until a few of years ago, all computations were performed in the quenched approximation 
in which the effect of virtual quark loops is neglected. Recent progress allowed us to move 
to the unquenched case in which at least Nf = 2 dynamical quarks are present in the QCD 
vacuum fluctuations. The lattice quark action that are being used nowadays are 0(a)-improved 
so that the systematic errors are O ((am c ) 2 ). In this write-up we present new results for the 
charmed decays using an improved Wilson action with Nf = 2 degenerate sea quarks. The 
results reported here refer to the simulations made at a ~ 0.08 fm for three different values 
of the sea quark mass corresponding to m n ~ 770 MeV, 600 MeV and 400 MeV.^ Other 
unquenched results relevant to the leptonic and semileptonic decays were obtained by using the 
staggered quark action^. Since the formal proof of validity of the staggered formulation is still 
missing, the study based on Wilson quarks is more than needed. 

2 Leptonic decays 

2.1 Hadronic matrix element 

The simplest way to determine the CKM matrix element \V c d\ is via the leptonic decay D — > iv 
with t = t, fi, e. The decay width is given by 

V(D+ -> Pv t ) = \v cd \ 2 ^m D+ mj (l - f D \ (1) 
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Figure 1: The extrapolation of the ratio fo/f-n either by using the HMxPT formula, eq. ([?]) (solid line) or a 

linear fit (dashed line). 

where Gp is the Fermi constant and the decay constant fn parametrizes the hadronic matrix 
element: 
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with = cj^j^d. The theoretical uncertainty in eq. ([I]) is entirely due to fjj. For I = e or [i 
this decay mode has been recently accurately measured^. On the lattice, /d is extracted from 
the asymptotic behavior of the 2-point Green function, i.e. 
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where is the appropriately renormalized axial current. Zq is evaluated from: 
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£ Computation of frj and chiral extrapolations 

We focus to the unquenched case in which the valence and the sea quark masses are equal. Thus, 
for each of our directly accessible light quark masses, we compute fu which then needs to be 
extrapolated to the physical fn d = Id- That extrapolation can be made either linearly in m q , 
or by using the expression derived in heavy meson chiral perturbation theory (HMxPT)^* 
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In this formula, $o an d c$ are the fit parameters, fo is the pion decay constant in the chiral 
limit determined on the same lattice, while g is related to the coupling between the heavy meson 
doublet, (D, D*) and the soft pion, i.e. 



(D(kMq)\D*(p,\)) = (e\q)g D . 
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were e A is the polarization vector of D*. Its value has been previously computed in the static 
limit - (g ~ 0.5), and with the propagating charm quark ^ (g ~ 0.6), both values leading 
to the large factor multiplying the logarithm in eq. ©, 3/4(1 + 3g 2 ). That would drive the 
extrapolation of fo q yJ m 'D q way below the result obtained through the linear extrapolation, thus 
to large systematic errors. It is then more convenient to consider the ratio /d/ fn in which the 
chiral logarithmic term is halved: 
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ci being a fit parameter. From fig. (pQ), we see that linear and x - l°g fits are very con- 
sistent and by using both values for g mentioned above, we get /d//tt = 1-50(24) j and 
Id/ fir = 1-52(17) linear . With the physical pion decay constant f n = 139.6 MeV^, we arrive at 

f D = 201(22) (±|) MeV . (8) 

3 Semileptonic decays 

3.1 Hadronic matrix element 

\Vcd\ can also be extracted by studying the partial or total decay width of the semileptonic decay 
D — ► irivi (i = e,n), which has been measured in various recent experiments^. The differential 
decay width is given by 

where q is the momentum transfer and X(q 2 ) = (q 2 + m 2 D — m 2 ) 2 — 4m 2 - ) m 2 . The vector form 
factor F + (g 2 ) parametrizes the hadronic matrix element of the weak current 

<k)\ (V - A)„ \D®) =(p+k- g W " F + (q 2 ) + q ^ ~ ^ F (g 2 ) , (10) 

where q = p — k and -F+(0) = -Fo(O). The contribution of the scalar form factor Fo (g 2 ) to 
the decay width comes with a mf -factor and therefore can be neglected. The form factors are 
extracted from the behavior of the following 3-point correlation functions: 

Cj?\k, q, t; Q = E 3 ,y (0| (™)o, (V»)y, t (576?)^ e<^) |o) 



computed on the lattice. The axial current does not contributed due to the parity conservation 

(3) 

in QCD. Extracting the matrix element from a fit of requires the knowledge of several 
independent quantities (-2^-,!); renormalization constants and 0{a) improvement coefficients of 
the vector current) whose error lower the accuracy of F+ (q 2 ). This can be improved by using 

the strategies we previously studied in^, based on ratios of such 3-point correlators and twisted 
boundary conditions. In fig. [2] we show the results obtained with the 1 st strategy. 

3.2 Computation of F + (g 2 = 1 GeV 2 ) and chiral extrapolations 

The differential decay width in eq. © is experimentally measured for various values of q 2 . To 
extract \ V c d\, one then needs a single F + (g 2 )-value. We choose q 2 = 1 GeV 2 where both theoret- 
ical and experimental errors are under a reasonable control. The extrapolation of F+(l GeV 2 ) 
toward the physical value leads to a difficulty similar to what we discussed above. To get 
around that difficulty, we use HM^PT fits for the ratio F + (g 2 ) / fjj where large deviations due 

to m 2 log m 2 terms are reduced^. We get 

GeV ) = 432 ( 56 ) GeV -i (HM X PT fit) ; 3.76(54) GeV" 1 (Linear fit) , 

ID 

where we also quote our result obtained by using the naive linear extrapolation. The difference 
of the two is an estimate of the systematic uncertainty of the extrapolation procedure. 
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Figure 2: The g 2 -dependence of the vector (left) and the scalar (right) form factors relevant to D — * n£v decay 
for 3 different pion masses, accessible directly from our lattices. 



4 Summary 

In this short note, we reported on the progress in determining the key hadronic quantities 
entering the leptonic and semileptonic decays on the lattice by using the Wilson quarks. A 
better control over the systematic uncertainties is achieved if one chooses judiciously the ratios 
in which various sources of uncertainties cancel out, or are diminished. In the case of semileptonic 
decays, also the use of twisted boundary conditions is very important. The quenched experience 
suggests that the 0(a 2 ) artifacts are reasonable^^, but to that end, more work is needed. 
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